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ABSTRACT 

It has been known for some time that the Milky Way is a barred disk galaxy. More recently several 
studies inferred from starcount observations that the Galaxy must contain a separate, new, flat long 
bar component, twisted relative to the barred bulge. Here we use a simulation with a boxy bulge and 
bar to suggest that these observations can be reproduced with a single structure. In this simulation a 
stellar bar evolved from the disk, and the boxy bulge originated from it through secular evolution and 
the buckling instability. We calculate starcount distributions for this model at different longitudes and 
latitudes, in a similar way as observers have done for resolved starcounts. Good agreement between the 
simulation and the observations can be achieved for a suitable snapshot, even though the simulation 
has a single boxy bulge and bar structure. In this model, part of the long bar signature is due to a 
volume effect in the starcounts, and another part is due to chosing a snapshot in which the planar 
part of the boxy bulge and bar has developed leading ends through interaction with the adjacent spiral 
arm heads. We also provide predictions from this model for the line-of-sight velocity distributions at 
the longitudes with the long bar signature, for comparison with upcoming surveys. 

Subject headings: Galaxy: structure — Galaxy: bulge — Galaxy: kinematics and dynamics 



1. INTRODUCTION 

During the last decades it has become clear that 
our Galaxy (hereafte r, the MW) is a barr ed system, 
as first suggested by Ide Vaucouleursl (119641). Now we 
have good evidence from NIR photometry dDwek et al. 
19951: iBinnev et al.l 119971) . star cou nts ()Stanek et al. 
1997b lLopez-Corredoira et al.l l2005[). g as kinematics 
(lEnglmaier fc Gerhardl 19991 : iFuxl I1999I ). microlensing 
fHamadach e et al.ll2006f) and dynamical effects near the 
solar circle ((Dehnen| 2000; Mi nchev et aTll2007l ). Several 
of these works refer to the boxy bulge which ends at a 
galactic radius of ~ 1.5 kpc , even though in so me models 
for the observed data (e.g.. IBinnev et aLlll997l ) it is clear 
that the MW bar extends further in the Galactic plane. 

Such a single boxy bulge and bar structure consist- 
ing of a boxy bulge and a planar bar continuation is a 
character istic outcome of the secular evolution of barred 
galaxies ()Athanassoulal 120051) . where th e bar eventu- 
ally buckles and forms a boxy-bulge (e.g.. iCombes et all 
119901: iRaha et al.l 119911 ). After the buckling event the 
bar resumes its evolution and continues to grow slowly 
through angular momentum exchange between bar, 
disk and dark matter hal o (e.g. iLvnden-Bell fe Kalnaia 
1 19721: lAthanassoulal I2003t iDebattista fc Sellwoodl l200of 
iMartinez- Valpuesta et al.ll2006l MV06). 

The MW's boxy bulge and bar has its long axis in the 
first quadrant, at an angle a ~ 15° — 30° wit h respect to 
the Sun-Galactic center line ([Ge rhard 2002). But more 
recently star count observations extendin g to greater lon- 
gitudes have led to surprising results flBcni ami n et al. 
2005 1 : ICabrera-Lavers et al.ll2007ll2008t iChurchwell et al. 
2009L here after B05,C07,C08,C09), c onfirming an earlier 
analysis bv lHammerslev et al.l (|2000| ): these observations 
found indications for a separate long bar with an orien- 
tation of a' ~ 43°, extending from longitude I ~ 27° to 
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I ~ 10° (4.5 kpc to 1.5 kpc galactocentric radius), and 
therefore coexisting with the conventional Galactic boxy 
bulge and bar over a range of radii. This interpretation, 
if correct, would dynamically be quite puzzling: two sep- 
arate rotating bars should align with each other through 
dynamical coupling in at most a few rotation periods. 

In this letter we show that a separate inferred long bar 
does not necessarily follow from the star count data, and 
we suggest a plausible model to explain these data with 
a single barred structure whose inner parts represent the 
boxy bulge. We also show some predictions for radial ve- 
locity distributions that could be used to test this model 
in the near future. 

2. A MODEL FOR THE MILKY WAY'S BAR AND BULGE 
FORMED THROUGH SECULAR EVOLUTION 

The simulation used in this work is similar to that pub- 
lished in MV06 and has not been run to match the MW 
struc ture. The code used is FT M 4.4 (updated version) 
from iHeller fc Shlosmanl (|1994l) . The total number of 
particles is 1 x 10 6 , distributed initially in an exponen- 
tial disk with Q = 1.5, embedded in a live dark matter 
halo. After ^1.5 Gyr the bar becomes very strong and 
buckles, thereby weakening. Later the bar resumes its 
evolution and grows again, resulting in a prominent boxy 
bulge and bar structure. 

We consider the simulated galaxy at time ~ 1.9 Gyr, 
after the boxy bulge has formed and the bar has regrown. 
The density distribution for this snapshot is shown in 
Figure Q]a, oriented at an angle a — 25° with respect to 
the line from the Galactic center to the observer. The 
boxy bulge is apparent in Figure [B). The model is scaled 
so that the end of the planar bar appears just inside 
longitude / = 30° as seen from the observer. The bar 
length is ~ 4.5 kpc, and the maximum ellipticity is 0.46. 
Relative to the planar bar, the boxy bulge is ~ 20% 
larger in I in the scaled model than in the MW, as mea- 
sured by comparing the maxima in the asymmetry map 
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Fig. 1. — Top panel: Face-on view of the simulation at time 
T ~ 1.9 Gyr. The bar rotates clockwise and its ends bend towards 
the leading side, connecting to the spiral arms further out. The 
model has been scaled to the MW and is oriented such that the 
long axis of the bulge is seen at an angle a = 25° by the observer 
at (0,-8 kpc). A second line at 43° as inferred for the long-bar in 
the MW is also shown. The two circles correspond to radii of 3 and 
4.5 kpc. The dotted lines show lines-of-sight for longitudes of 0°, 
10° , 20° and 30° from the observer. Lower panel: edge-on view of 
the same snapshot, as viewed from the Sun. The boxy structure is 
noticeable. Higher densities correspond to brighter colors. 

with those from the COR E data (jBissantz et al.l 119971; 
IBissantz fc Gerhard! I2002D . Also both the boxy bulge 
and the disk are vertically more extended in the model. 
The Sun is placed at 8 kpc. 

In the face-on view we can easily identify the curved, 
leading ends of the stellar bar. Over a period of 
1.2 Gyr, the model shows oscillations from leading 
through straight to trailing ends and back. The bar 
spends 40% of this time in the leading phases. Simi- 
lar morphology can be seen in o ther barred simulations 
in the literature (e.g. iFuxl 119971 model m08) and also 
in so me observed gal axies such as NG C 3124 (|Efremovl 
l20ll and NGC 3450 (jButa et al.ll2007h . The oscillations 
between trailing and leading ends of the bar could be re- 
lated to the oscill ations seen in the bar growth in N-body 
simulations (e.g. IDubinski et al J 120091 ) and may be due 
to no n-linear coupling m odes between the bar and spiral 
arms (jTagger et al.lll987ft . This topic is beyond the scope 
of this paper. For comparison we show a snapshot at a 
later time in this simulation where the ends of the bar 
are straight and the spiral arms appear to emerge from 
them (Figured]). 
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Fig. 2. — A snapshot from the same simulation, but at this par- 
ticular time the ends of the bar are aligned with the bar's inner 
part. This snapshot is used to illustrate the volume effect but is 
not as good a representation of the MW as that in Figure 1; see 
Section 2.1. 

2.1. Quantitative analysis: no separate 'long bar' is 

needed 

We apply a similar technique as was used to iden- 
tify two different barred structures in the MW from star 
count data (B05,C07,C08,C09). We view the projected 
model as an observer in the disk at 8 kpc distance from 
the center would see it. To increase the particle resolu- 
tion, we symmetrize the model vertically and divide the 
(longitude, latitude)-{l , &)-space into bins of 51 = 3° and 
5b = 2°, respectively. Then we count particles in each 
of the corresponding cones and bin these particles in dis- 
tance modulus, with S[i = 0.1. The distance modulus is 
given by [i = —5. + 5. x log(D[pc\). 

In Figure [3] we show histograms in n for several lines- 
of-sight. To quantify the distribution of particles with n 
and to assign a distance value to the maximum number 
counts, we fit a Gaussian to the left- most peak, i.e., the 
one nearest to the observer. In the histogram obtained 
when looking towards the ends of the bar in the Galac- 
tic plane (Figure \5ja) we can identify three main peaks, 
one corresponding to the bar, one to a spiral arm in the 
back, and one to the end of the disk. The /z-value of the 
fitted maximum corresponds to the end of the bar, where 
the bar is flat. The second histogram for (I, b) — (9°, 8°) 
(Figure shows the distribution of stars in a field well 
above the plane where the boxy bulge dominates. The 
fitted maximum corresponds to a position on the thick 
line in Figure [T] at a = 25°. In the first panel we can 
clearly identify the particles in the disk, but in the sec- 
ond showing the higher latitude field, disk particles are 
absent. We also show the histogram for (9°,0°) (Fig- 
ure^), where we can see the increment in the number 
of particles with respect to those at (27°, 2°), and the 
displacement of the maximum towards larger distance in 
comparison with (9°, 8°). In the last panel we show one 
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Fig. 3. — Distribution of particles with distance modulus in four 
fields as seen by an observer at the Sun's position 8 kpc from 
the center. The left panels a), c) show distributions in/near the 
Galactic plane in cones centered on the given longitude and lati- 
tude. The right panels b), d) show histograms for cones through 
the boxy bulge. 

of the central lines-of-sight (3°, 6°) for the boxy bulge. 

Repeating the identification of maxima for all fields in 
-6° < I < 30° and -8° < b < 8°, we produce a similar 
plot as in C07. The distance modulus for each maxi- 
mum is converted to distance in kpc from the observer's 
position and plotted in Figured] in galactocentric coor- 
dinates. We show such a plot both for the snapshot with 
leading bar ends (Figure [Da, corresponding to Figured]), 
and for the snapshot with straight ends (Figure |4^, corre- 
sponding to Figure [2]). In both models, the points com- 
puted from star counts above the plane (pink crosses) 
follow the heavy line indicating the true orientation an- 
gle a = 25° fixed by us for the model's boxy bulge and 
bar, except at the very center and at negative I where the 
maxima are closer to the observer due to tangent point 
effects. 

For lower latitudes (black crosses), the points move fur- 
ther back. In both models, for longitudes 0° < I < 9° in 
the region of the boxy bulge, points in the plane now ap- 
proximately follow an imaginary line with a' = 43°. This 
is due to a volume effect combined with the shallow disk 
density distribution: N(n) oc n{D)D 2 dD / 'd/i. There is 
an intermediate region with 9° < I < 20°, where the 
boxy bulge becomes thinner and transits into the planar 
bar, and then there is a region, 20° < I < 30° (~ 3 to 
~ 5 kpc), where we can see just particles in the plane. In 
this region, both models differ: the volume effect alone, 
as illustrated for the bar with the straight ends in Fig- 
ure [4]!), accounts for only half the shift from a = 25° to 
a' = 43°. Whereas for the snapshot with the leading bar 
ends (Figure |4]a) the in-plane points now approximately 
follow the imaginary a' = 43° line right until the end 
of the bar is reached. If we were not aware of the origi- 
nal structure of the model, the results shown in Figure |4] 
could easily be interpreted as evidence for a model with 
two structures, one thicker and shorter at a — 25° and 
another thinner and longer at a' ~ 43°. Clearly, the case 
discussed shows that a separate long bar is not necessar- 
ily implied by the observed star count distributions and 
that a single boxy bulge and bar structure is a valid and 
dynamically simpler interpretation of these data. 
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Fig. 4. — Location of the star count maxima in the Galactic 
plane, for fields near the disk plane (black crosses) and in the boxy 
bulge (4° < \b\ < 8°, pink crosses). The top panel (a) shows the 
maxima for the model with leading curved ends of the bar; black 
dots show the maxima for the initial axisymmetric disk. The lower 
panel (b) is for the model with straight bar ends. The circled 
crosses correspond to the histograms shown in Figure \3\ The thick 
solid line shows the true orientation of the model, a = 25°. The 
thin line follows a' = 43° . The dashed circle has a radius of 3 kpc. 

3. RADIAL VELOCITY PREDICTIONS 

Several large radial velocity surveys of Galactic bulge 
and disk stars are currently on-going or planned, and 
we thus provide some predictions for our model. We 
scale the a verage cir c ular v elocity curve of the model to 
the data of lClemensI (|1985() for radii 3-4.5 kpc, and then 
compute the radial velocities at different lo ngitudes as 
seen from the Sun (assuming w Q =250 km/s: iReid et al.1 
2009]). We show the mean velocity, dispersion, and par- 
ticle number vs. distance modulus in Figure [S] for two 
in-plane positions near the end of the planar bar. The 
non-circular motions in the barred models leave a clear 
signature at ji ~ 13.0 — 13.6 compared with the initial 
axisymmetric disk. In this range the bar's velocity pro- 
file is fiat. For the leading-ends bar it has a break at 
/i ~ 13.7 and then increases towards an axisymmetric 
type of profile. For the straight-ends bar the break hap- 
pens earlier /i ~ 13.5 and the radial velocities are slightly 
smaller. The difference between axisymmetric and bar 
case is ~ 30 — 40 km/s. The velocity dispersions of both 
bar models in the two fields are higher by ~ 20 km/s for 
all /U < 14.2. Also note the difference between the star 
counts with ^ between the three models. 

Velocities in the boxy bulge are outside the scope of 
this paper; they have already been studied with a similar 
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Fig. 5. — Dependence of mean radial velocity (top), velocity dis- 
persion (middle), and particle number (bottom) on distance mod- 
ulus, for two fields centered at 27° and at 24°. Solid, dotted, and 
dot-dashed lines are for the snapshot with leading-ends of the bar, 
for the bar with straight ends, and for the initial exponential ro- 
tating disk, respectively. The signature of the bar can be easily 
identified; see text. 

N-body simulation by IShen et all (|2010D and compared 
with results from the BRAVA survey. 

4. SUMMARY AND CONCLUSIONS 

We have analyzed star counts in the inner Galaxy using 
an N-body model which arose from secular evolution of 



a disk galaxy, a natural mechanism for galaxies like ours. 
The bar and buckling instabilities in the stellar disk lead 
to a boxy bulge which extends to a longer in-plane bar. 
The bar couples with the spiral arms in the disk, giving 
rise alternately to leading, straight or trailing bar ends. 

As seen from within the disk at 8 kpc from the cen- 
ter, the maxima of the line-of-sight distance distributions 
in the Galactic plane occur at distances somewhat fur- 
ther than the maxima of the line-of-sight density dis- 
tributions, due to the volume effect in the star counts. 
Assuming a plausible orientation (a = 25°), this explains 
part of the observational signature which was previously 
used to infer the existence of a second long bar. If in 
addition we choose a model snapshot where the bar has 
leading ends, most of the long bar signature in the star 
count data can be reproduced. While not made specially 
to match the MW, this model thus illustrates that the 
traditional Galactic bar (the boxy bulge) and the more 
recently inferred long bar can plausibly be explained by 
a single boxy bulge and bar structure. 

To test this further we have determined the depen- 
dence of the mean radial velocity and velocity dispersion 
on distance modulus in Galactic plane fields near the 
inferred end of the planar bar. These illustrate the dif- 
ferences between the barred model and an axisymmetric 
rotating disk, which can be compared with upcoming ra- 
dial velocity survey data for the inner MW. Future work 
should also address in more detail the spiral arm - bar 
interaction which gives rise to the curved ends of the 
bar, and aim at constructing a more detailed dynamical 
model allowing us to understand better the structure and 
evolution of the inner Milky Way. 
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